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Abstract 

We study spherically and axially symmetric monopoles of the SU(2) Einstein- 
Yang-Mills-Higgs-dilaton (EYMHD) system with a new coupling between the 
dilaton field and the covariant derivative of the Higgs field. This coupling arises 
in the study of (4 + 1)- dimensional vortices in the Einstein- Yang- Mills (EYM) 
system. 
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1 Introduction 



In the 1920s, Kaluza and Klein studied the 5- dimensional version of the Einstein 
equations Q by introducing a 5- dimensional metric tensor. When one dimension is 
compactified, the equations of 4-dimensional Einstein gravity plus Maxwell's equa- 
tions are recovered. One of the new fields appearing in this model is the dilaton, a 
scalar companion of the metric tensor. In an analog way, this field arises in the low 
energy effective action of superstring theories and is associated with the classical scale 
invariance of these models 0. 

In recent years, a number of classical field theory models coupled to a dilaton have 
been studied. It was found that if the solutions of SU(2) Yang-Mills-Higgs (YMH) 
theory, namely the 't Hooft-Polyakov monopole H and its higher winding number 
generalisations @, [|, U, are coupled to a massless dilaton [|7|, remarkable similarities 
to the qualitative features of Einstein- Yang-Mills-Higgs (EYMH) monopoles |9|, [10 



arise. Especially, it was observed that similarly to gravity the dilaton can render 
attraction between like charged monopoles and thus bound multimonopole states are 
possible. 

Recently, the YMH system coupled to both gravity and the dilaton has been stud- 
ied g. It was found that in the full Einstein- Yang-Mills-Higgs-dilaton (EYMHD) 
model, a simple relation between the tt-component of the metric and the dilaton ex- 
ists. The abelian solutions to which the configurations tend in the limit of critical 
coupling are the extremal Einstein-Maxwell-dilaton (EMD) solutions. These have 
definite expressions for the energy and the value of the rr-component of the metric 
at the origin which only depend on the fundamental couplings. 

Volkov argued recently [BJ that if is a symmetry of the Einstein- Yang-Mills 
(EYM) system in (4 + 1) dimensions, where £4 is the coordinate associated with the 
5th dimensions, then the (4 + 1) -dimensional EYM system reduces effectively to a 
(3 + l)-dimensional EYMHD system with a specific coupling between the dilaton field 
and the Higgs field. 

In this paper, we study both the spherically and axially symmetric solutions of 
the (3 + l)-dimensional EYMHD model deduced from the (4 + l)-dimensional EYM 
system. In Section 2, we give the EYMHD Lagrangian and review how the specific 
coupling arises. In Sections 3 und 4, we give the Ansatz and present our numerical 
results for the spherically and axially symmetric solutions, respectively. In both these 
Sections emphasis is placed on the flat space limit and on the limit in which our system 
of equations reduces to the one in [jEJ. The summary and conclusions are presented 
in Section 5. 



2 SU(2) Einstein- Yang-Mills-Higgs-dilaton theory 

The Lagrangian and the particular coupling of the dilaton field \& to the SU(2) gauge 
fields and Higgs fields $ a (a = 1,2,3), respectively, arise effectively from the 
Lagrangian of (4 + l)-dimensional EYM theory. If both the matter functions and the 
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metric functions are independent on x 4 , the 5-dimensional fields can be parametrized 
as follows (with M, N = 0, 1, 2, 3, 4) Q: 

g%dx M dx N = e-tgWdx^da? - e 2( {dx 4 ) 2 , n, v = 0, 1, 2, 3 (1) 

and 

A^x M = + $ a da; 4 , a = 1, 2, 3 (2) 

where is the 4-dimensional metric tensor and ( plays the role of the dilaton. 
Introducing a new coupling k to study the influence of the dilaton systematically, we 
set ( = 2k\I / and obtain the following action of the effective 4-dimensional EYMHD 
theory: 

S = S G + S M = J L G ^g~^)d A x + J L M ^g~^)d A x . (3) 
The gravity Lagrangian Lq is given by 

iG = itd R • (4) 

where G is Newton's constant, while the matter Lagrangian Lm reads: 

L M = F^F^ - -d^d^ - -e-^D^ a D^ a - e- 2 ^V{^ a ) , (5) 

1 2 2 



with Higgs potential 
the non-abelian field strength tensor 



V($ a ) = -(<$> a ® a -v 2 ) 2 , (6) 



F; v = d,A a u - d v Al + ee ahc A\Al , (7) 

and the covariant derivative of the Higgs field in the adjoint representation 

D^ a = 8^ a + ee abc A b ^ c . (8) 

Here, e denotes the gauge field coupling constant, A the Higgs field coupling constant 
and v the vacuum expectation value of the Higgs field. 



3 Spherically symmetric solutions 

For the metric, the spherically symmetric Ansatz in Schwarzschild-like coordinates 
reads |§: 

ds 2 = g$dx»dx v = -A 2 {r)N{r)dt 2 + N~\r)dr 2 + r 2 d6 2 + r 2 sin 2 0d?<p , (9) 

with 

N(r) = 1-^1. (10) 
r 
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In these coordinates, m(oo) denotes the (dimensionful) mass of the field configuration. 
For the gauge and Higgs fields, we use the purely magnetic hedgehog ansatz [|3| 

AS = A t a = , (11) 

A a = ^~K{r) = _ 1-K(r) ^ ^ a 

e e 

$ a = v H{r)e r a . (13) 

The dilaton is a scalar field depending only on r 

^ = $(r) . (14) 

Inserting the Ansatz into the Lagrangian and varying with respect to the matter fields 
yields the Euler-Lagrange equations, while variation with respect to the metric yields 
the Einstein equations. 

With the introduction of dimensionless coordinates and fields 

$ m 

x = evr , u = evm , (/) = — , ip = — . (15) 

v v 

the Lagrangian and the resulting set of differential equations depend only on three 
dimensionless coupling constants, a, (3 and 7, 

r- M w „ V\ M H kM w 

a = vGv = — — — , a = — = — = , 7 = kv = , 16 

eM P1 H e ^/2M W e V ; 



where M w = ev, M H = V2Xv and M Pi = 1/VG. With the rescalings (0) and flT6|), 
the dimensionless mass of the solution is given by 

With ( jl5|) and ([U]) the Euler-Lagrange equations read: 

(e 2 ^ANK'Y = A(e 2 ^ K ^ K \~ l) + e~^H 2 K) , (17) 

x 2 

(e-^x 2 ANH')' = AH(2e~ 4 ^K 2 + f3 2 x 2 e~ 2 ^(H 2 - 1)) , (18) 

(x 2 AN^J = 2 1 A[e 2 ^(N(K') 2 + ^— ^) 

2x 2 

_ e -2^f^( H 2 _ 1)2 _ 2e -^{-N{H') 2 x 2 + H 2 K 2 )} , (19) 

where the prime denotes the derivative with respect to x, while we use the following 
combination of the Einstein equations 

G tt = 2a 2 T tt = -2a 2 A 2 NL M , (20) 
9 XX G XX - g*G tt = -4« 2 iV^ (21) 
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to obtain two differential equations for the two metric functions: 



// = a 2 (e 2 ^N(K') 2 + -Nx 2 {H') 2 e~^ + — 2 {K 2 - 1) + K 2 H 2 e~ A ^ 

+ ^x 2 (H 2 - 1)V 2 ^ + ^iVxV) 2 ) , (22) 

A' = a 2 xA (^^Le 2 ^ + e-^{H') 2 + (^) 2 ) . (23) 

Since we are looking for globally regular, finite energy solutions which are asymptot- 
ically flat, we impose the following set of boundary conditions: 

tf(0) = l, #(0) = 0, d x ij>\ x=o = 0, /i(0) = 0. (24) 

fsT(oo) = , H(oo) = 1 , V(oo) = > A (°°) = 1 > (25) 



3.1 Numerical results 

We have restricted our numerical calculations for the spherically as well as for the 
axially symmetric solutions to (3 = 0. 



3.1.1 The a = limit 

For a = 0, the gravitational field equations (|22"D and (p3|) decouple from the rest of the 
system and we are left with the YMHD model in flat space, N(x) = 1 and = 1. 
In it was observed that (in analogy to the EYMH system) the monopoles exist up to 
a maximal value of the dilaton coupling 7 = 7 ma x and from there on a second branch 
of solutions tend to the abelian solution for 7 — * 7 cr < ^max with ^(0) monotonically 
decreasing to —00. Our numerical results indicate that no such 7 ma:r exists in the 
model studied here. We have integrated the equations for 7 e [0 : 10] and found the 
solutions to exist for all these values of 7. The profiles of the functions suggest that 
for 7 — > 00 the solutions tend to the vacuum solution with K(x) = 1, H(x) = 
and ijj(x) = 0. This is indicated by the fact that the mass of the field configurations 
progressively tends to zero for rising 7 and that the matter fields are equal to their 
vacuum values on increasing intervals of the coordinate x. To demonstrate this, we 
give below the values of xk and xh, where the gauge field K(x) and the Higgs field 
H(x), respectively, reach the value 0.5, i.e. K(xk) = 0.5 and H(xh) = 0.5 : 



7 


x K 


x H 


0.1 


2.2 


1.8 


1.0 


3.3 


2.1 


5.0 


11.2 


5.1 


10.0 


22.4 


8.9 



5 



Moreover, we find that ip(0) > for all 7. Since for 7 = 0, the BPS monopole solution 
is recovered, the curve for ^(O) starts from zero at 7 = 0. From there it increases to 
a maximal value at 7 ~ 1.38 and then slowly decreases to zero for 7 — > 00. Together 
with (p5|) our results strongly suggest that ip(x) tends to zero on the full interval 
x e [0 : 00 [. 

The difference between the model studied here and the one studied in is that in 
the standard case of the Yang-Mills-Higgs-dilaton system, the equations are effectively 
the equations of an Einstein- Yang-Mills-Higgs model with metric 

ds 2 = -e 2 ^dt 2 + e- 2K ^dr 2 + e - 2K ^r 2 d6 2 + e~ 2 ^\ 2 sin 2 9d 2 V (26) 

The components of the Einstein tensor then read: 

Goo G 



G tt = -Ke 



_ 4K*(r) 



2$"+ 4— -k(&) 2 
r 



■ Grr = = = «W (27) 

r z r z sm 6 



The following combination of the Einstein equations 

G\ - G r r - G e e - G% = MGiTl - T; - T e e - T«) (28) 

exactly gives the dilaton equation in the standard case with 7 2 = k 2 v 2 oc Gv 2 . 
Equally, the Euler-Lagrange equations for the gauge and Higgs field functions are 
obtained using the metric (p6|) . Thus the dilaton can be viewed as a metric field. 
Since we know that gravity leads to a critical value of the coupling constant, this is 
also true in the standard Yang-Mills-Higgs-dilaton case. A horizon forms for g tt — ► 
which implies \& — > —00. This is exactly the limiting solution found previously. 

Now, for the equations studied here, this is not possible. We cannot introduce a 
4-dimensional metric which makes the equations in the limit of A = N = 1 reduce 
to a set of Einstein- Yang-Mills-equations. yThus, no critial 7 can be expected on the 
basis of the above argument. This is confirmed by our numerical results. 

The form of the Lagrangian (|5|) suggests that for 7 getting bigger and bigger, the 
only possibility to fulfill the requirement of finite energy is that ip{x) — > on the full 
interval of x. 



3.1.2 The Volkov limit a 2 = 3-f 2 

In |L2| only one fundamental coupling is given, the gravitational coupling G. Com- 
paring the (4 + l)-dimensional EYM system with the EYMHD system studied in this 
paper, we conclude that for 

a 2 = 37 2 (29) 

our system of equations reduces to the one in [0 . 

Studying the full system of equations, we were particularly interested in reobtain- 
ing the results with a different numerical method |13| and in studying some of the 
features of the solutions in greater detail. 
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Volkov observed a spiraling behaviour of the parameters for the EYM vortices. As 
expected and shown in Fig. 1, we observe this feature as well. Our numerical results 
suggest that a number of branches exist, on which the minimum of the metric function 
N, N m , and the value of the metric function A at the origin, A(0), monotonically 
decrease. In the limit of critical coupling, the rr-component of the corresponding 
5-dimensional metric tensor, expressed in 4 dimensions through the function g(x) : = 

e 2 ^N(x) = e~J3 a ^ N(x) develops a double zero at a value x m > of the dimensionless 
coordinate x |T2|]. Our numerical results indicate that x m e ]0 : 0.082]. 

There exist several local maximal and minimal values of a, ac m ^ x and oc^ n , respec- 
tively. By a^xmin we denote the value of a, at which the i-th and the j-th branch 
join. We find: 

ofi2, = 1.268, offi = 0.312, ag£ = 0.419, «g£ = 0.395 (30) 

Apparantly, the difference between a max and its corresponding a min decreases, i.e. 
the branches get smaller and smaller. We thus conjecture that the limiting solution 
is reached for a value of a close to ot^ n = 0.395. 

Progressing on the branches, the qualitative behaviour of the functions changes 
I"2fl . This is shown for the gauge field function K(x) in Fig. 2. While for all values 



of a on the first branch (as an example, we show K(x) for a = 1.267 close to cx^x) 
and most of the values of a on the second branch (see the profile for a = 0.8) K{x) 
decreases monotonically from its value at the origin K(0) = 1, to its value at infinity 
K(oo) = 0, oscillations of the gauge field start to occur on the second branch for 

(23) 

a < 0.4185. For or mi ^ = 0.312 the local minimum and maximum, respectively, are 
already quite pronounced, while for a^f„ = 0.395, the location of both the minimum 
and maximum has moved to smaller values of x. We are convinced that in analogy 



to [12| the number of oscillations increases when proceeding on further branches. In 
Fig. 3, we show the value of the dilaton function if> at the origin (multiplied by — 1), 
— ^(O), and the mass of the solutions on the different branches as functions of a. 
ip(0) is equal to zero for a = (since this implies 7 = 0), increases to a maximal 
value of ^(0) = 0.174 at a ~ 1.0 and from there decreases first to zero at a = 1.220 
on the second branch and further decreases to a minimal value of ip(0) = —5.274 at 
a = 0.314 on the third branch. From there it starts to increase to ip(0) = —4.847 at 
a = 0.395 and then decreases again to ^(0) = —7.164 at = 0.395. 

The function ip(x) itself decreases monotonically from ip(0) to its value ^(oo) = 
for all values of 7 on the first branch and for a < 1.24 on the second branch, staying 
positive for all values of x. For a > 1.24 on the second branch, a minimum (which 
has negative value) starts to form and dips down deeper to negative values when 
progressing on the branches. In the limit of critical coupling a cr , the dilaton function 



i>EMD of the corresponding EMD solution (TT[ 



^MD = ^-ln(l-^), X_ = (V* (31) 
Aa rr X 3 
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where the dimensionless coordinate x is given in terms of X , 



x = a cr X(l - (32) 

is reached for x > x m , while it stays finite for x e [0 : x m [. 

It has to be remarked here, though, that in the limit of critical coupling, the gauge 
and Higgs field functions K(x) and H(x) don't reach there abelian values and 1 at 
x = x m , respectively, but tend to the fixed point described in |T2"]. 



The mass of the solution stays close to one on the first branch and increases 
monotonically on the second branch. On the third and fourth branch it differs only 
little from the corresponding mass on the second branch and thus the three different 
branches can barely be distinguished in Fig. 3. In the limit of critical coupling a cr , 
the mass tends to the mass ^ Et f D of the corresponding EMD solution: 



H _^ Hemd 




(33) 



Let us finally investigate the Einstein- Yang-Mills-Higgs (EYMH) system with the 
usual coupling of the dilaton considered in JTT| in the Volkov limit a 2 = 3J 2 . This 



model doesn't contain the prefactor of the covariant derivative involving the dilaton 
studied here. In the BPS limit (/? = 0), the equation for the Higgs field doesn't 
involve the dilaton field directly and vice versa, while in the model studied here it 
does even for (3 = (see ( |18D and (|I~9D ). In Fig. 1, we show the values of N m and 
A(0) for the model studied in [|11[] for the Volkov limit a 2 = 37 2 . The recalculation 
of the configurations for this specific relation between a and 7 exactly confirms the 
results obtained in |rTJ. The solutions exist for a < a max = 1.216, which fulfills the 



condition (48) of fTl 



«L + 7 2 = y g<W « 1.4 . (34) 

From there, on a second branch of solutions, the Einstein-Maxwell-dilaton (EMD) 
solution is reached at a cr with the minimum of the metric function N, N m , tending 



to the value N(x = 0) of the corresponding EMD solution ||15|| : 



N m -> N EMD (0) = i-^-^) 2 = 4 (35) 
a z + 7 2 lb 

This is demonstrated in Fig. 1. A(0) decreases to zero, while N m > 0, indicating that 
the extremal EMD solutions have a singularity at x = not hidden by an horizon. 



4 Axially symmetric solutions 



Since the n = 1 monopole was shown to be the unique spherically symmetric solution 
in SU(2) YMH theory [TBI], we need to impose an axially symmetric Ansatz (or one 



S 



with even less symmetry) to construct higher winding number solutions. The axially 
symmetric Ansatz for the metric in isotropic coordinates reads: 

ds i = - fd e + ^ (dr 2 + f 2 d9 2 ) + j? 2 sin 2 9dip 2 , (36) 

The functions /, m and I now depend on f and 9. If Z = m and / only depend on f, 
this metric reduces to the spherically symmetric metric in isotropic coordinates and 
comparison with the metric in (|9]) yields the coordinate transformation 

dr 1 dr 



r sjN(r)r 

For the gauge fields we choose the purely magnetic Ansatz ||: 



(37) 



A t a = , A f a = ^ < , (38) 
er 

A e a = l —^v/ , V = -- sin 9 ( H ^r a + (1 - H 4 )v e a ) . (39) 

e e 

while for the Higgs field, the Ansatz reads [[5], 

= v (<$> lVf a + $ 2 ^ a ) . (40) 

The vectors v?,vo and are given by: 

Vf = (sin9 cos nip, sin 9 sin nip, cos 9) , 
vg = (cos 9 cos nip, cos 9 sin nip, — sin 9) , 

= (— sin nip, cos nip, 0) . (41) 

The dilaton field \I/ now depends on f and [Q: 

# = tf(f,0) (42) 

For H 1 = H 3 = ^ 2 = 0, H 2 = H A = K(r), $x = /7(f), * = *(f) and n = 1, this 
Ansatz reduces to the spherically symmetric Ansatz in isotropic coordinates. 

The Euler-Larange equations arise by varying the Lagrangian with respect to the 



matter fields, while we use the following combinations of the Einstein equations [19 



grf{G^-2G t ) = 167rG T {L M + g — —-g -—r - g ——- g ^ ) (43) 



f y * dg u * dg ff * dg ee * dgw 
g ff {Gt + Gl) = l6nG^(L M - g^ - g^) (44) 

g fF (Gt + G 6 ) = 1GkGJ(L m - g^ - g^) (45) 



to obtain 3 differential equations for the metric functions /, I and m, which are 
diagonal with respect to the set of derivatives (/,f,f, ??V,f, hr,f-> l,9,0i With an 

analog rescaling as in (|i5|), the set of 10 partial differential equations again depends 
only on the three fundamental coupling constants introduced in ([16]). 
At the origin, the boundary conditions read (with 8% = -^d?)\ 

d s f(0,9) = d i l(p,9)=d £ m(p,9) = 0, akV(O,0) = O (46) 

Hi(0, 0) = O,i=l, 3, Hi(0, 0) = l,i = 2, 4, 0,(0, 9) = 0, % = 1, 2 (47) 

At infinity, the requirement for finite energy and asymptotically flat solutions leads 
to the boundary conditions: 

/(oo, 9) = l(oo, 9) = m(oo, 9) = 1, ^(oo,^) = (48) 

Hi(oo, 9) = 0, i = 1, 2, 3, 4, ^(oo, 0) = 1, 2 (oo, 0) = (49) 

In addition, boundary conditions on the symmetry axes (the p- and z-axes) have to 
be fulfilled. On both axes: 

H x = H 3 = 2 = (50) 

and 

d e f = d e m = d e l = d e H 2 = d e H 4 = d e <pi = d e ip = (51) 
4.1 Numerical results 

Constructing the axially symmetric solutions numerically [P0| , we were mainly inter- 



ested whether bound multimonopoles are possible in this model and if so, how the 
strength of attraction compares to that in ||. 

4.1.1 The a = limit 

Since gravity itself is known to be attractive, we first studied the influence of the 
dilaton alone. Constructing the n = 2 solutions, we obtain the same qualitative 
feature as for the n = 1 case. The configurations exist for all values of 7, tending 
to the vacuum solution Hi = H% = <pi = <p2 = ip = 0, H2 = H4 = 1 for 7 — ► 00. 
Again, we find that ip(0) > and that the curve for "0(0), starting from zero at 7 = 0, 
develops a maximum at 7 = 1.38 (which coincides with the corresponding value for 
n = 1). From there it monotonically decreases to zero for 7 — > 00. 

In Fig. 4, we show the difference between the mass of the n = 1 solution and the 
mass per winding number of the n = 2 solution AE = E(n = 1) — E(n = 2)/2. The 
sign of AE indicates whether like charged monopoles attract or repell each other, 
the modulus of AE indicates the strength of attraction and repulsion, respectively. 
We find that for all values of 7 we have studied, the monopoles are in an attractive 
phase and that with increasing 7 the strength of attraction grows. For the model 
studied in ||, the maximal AE was found to be AE max w 0.0082 at 7 = 1.2. From 
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there AE decreases since the spherically as well as the axially symmetric monopoles 
are tending to an essentially abelian, spherically symmetric solution in the limit of 
critical coupling. In the model studied here, AE is smaller for the same values of 7 
(e.g. we find that AE{j — 1.2) pa 0.0069), but since no critical coupling exists, AE 
grows with increasing 7. For e.g. 7 = 5.0, we find AE pa 0.0167, which is more that 
twice the maximal possible value for the model in ||. For large enough values of 7, 
AE should decrease, because both the spherically and axially symmetric solution are 
tending to the vacuum. Since the absolute value of the mass is decreasing to zero, 
the same should happen for AE. 



4.1.2 The Volkov limit a 2 = 3-f 2 

Studying the n = 2 axially symmetric monopoles, we find the same qualitative fea- 
tures as in the spherically symmetric case. We obtain: 

affi = 1.279, a$l = 0.295 (52) 

These values are slightly bigger and smaller, respectively, than the corresponding 
values for the n = 1 solution. It has already been observed previously, that the 
branches for the n = 2 solutions exist for higher values of the maximal coupling [IT? 



We find here that an analog holds true for the minimal coupling, namely that the 
n = 2 solutions exist for smaller minimal values than in the n = 1 case. 

The behaviour of ^(0) is completely analog to the n = 1 case : starting from zero 
at a = 0, it increases to a maximal value of ip(0) ~ 0.178 at a ~ 1.0, then decreases to 
zero at a ~ 1.19 and reaches ip(0) pa —5.324 at or^ n = 0.295. The dilaton function, 
which has only a weak angle-dependence, is positive and monotonically decreasing on 
the first branch and starts to form a negative valued minimum on the second branch. 
Proceeding on the second branch, the solutions are less and less angle dependent 
and the maximum of the function Hi, if 3 and $ 2 decreases. We didn't manage to 
construct further branches with our numerical routine (since these branches are rather 
small), but believe that in analogy to the n = 1 solutions, the n = 2 solutions develop 
a double zero, which is now (because of the choice of isotropic coordinates) located 
at x = 0. 

Comparing AE for the model studied here and for that in ||, we find that the 
different values don't deviate very much from each other. However, in the model 
studied in |§ the parameter space, in which solutions exist, is limited by \/a 2 + 7 2 ~ 
1.4. This means that in the Volkov limit solutions only exist for a < a max ~ 1.21. 
Here, the solutions exist for bigger values of a = v^3~7 and thus the maximum of AE 
is slightly bigger with AE max pa 0.011. This is reached on the first branch of solutions 
as demonstrated in Fig. 4. The second branch starts from a value of AE well below 
the corresponding one on the first branch and from there decreases. For all values of 
7 = the value of AE on the first branch exceeds the value of AE on the second 
branch. 
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5 Conclusion and Summary 



We have studied spherically and axially symmetric solutions of SU(2) EYMHD theory. 
The specific coupling between the dilaton and the Higgs field in this model arises from 
a (4 + 1) -dimensional EYM model studied recently in [12]. We find that in the flat 



space limit (a = 0), the solutions exist for all values of 7 tending to the vacuum 
solution for 7 — > 00. This contrasts to the situation in the usual YMHD model 
studied in 0, where the solutions only exist for 7 < 7 max tending to an abelian 
solution in the limit of critical coupling 7 cr . In both models, the dilatonic monopoles 
are in an attractive phase. Since no critical coupling exists in the model studied here, 
much stronger bound monopoles are possible. 

Apparently, unlike the YMHD model with the usual coupling, the YMHD system 
studied here doesn't have analog features than the EYMH system. Moreover, the 



relation between the tt-component of the metric and the dilaton found in ||11|| , doesn't 
exist 



When in addition gravity is coupled, our equations agree with those in |12j for 
a specific relation between the gravitational coupling a and the dilaton coupling 7, 
namely a 2 = 37 2 . Studying the spherically and axially symmetric solutions of these 
equations in this limit, we observe the same spiraling behaviour of the parameters 



found in [ 12j] . Several branches exist differing from the situation in the model with 



the usual dilaton coupling studied in [11]], where maximally two branches exist. Like 



in [H the dilatonic monopoles can form bound multimonopole states. 



In [|TTJ, non-abelian black holes of the EYMHD model were studied giving further 
evidence that the "No-hair" conjecture doesn't hold in models involving non-abelian 
fields. It would be interesting to construct the analogs in this system, especially the 
axially symmetric non-abelian black holes. Their existence would show that Israel's 
theorem cannot be extended to EYMHD theory, like it has previously been observed 
in EYMH theory 0. 

Finally, since in the low energy effective action of string theory, further correc- 
tions to the equations of standard physics, such as higher order corrections to gravity 
in form of Gauss-Bonnet terms or new fields like the axion field arise, it would be 
interesting to study their influence on the solutions constructed here. 

Acknowledgements B. H. was supported by the EPSRC. 
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Figure 2: The gauge field function K(x) is shown as function of the dimensionless 
coordinate x for the n = 1 solutions in the Volkov limit a 2 = 3"y 2 for the following five 
values of a : 1 : a = 1.267 close to oc^ax , 2 : a = 0.8 (second branch), 3 : a = 0.4185 
(second branch), 4 : a^l = 0.312 and 5 : a^ n = 0.395. 
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Figure 3: The value of the dilaton function tp at the origin (multiplied by —1), -t/>(0), 
is shown as function of a for the n = 1 solutions in the Volkov limit a 2 = 37 2 . Also 
shown is the mass of these solutions. 
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Figure 4: The difference between the mass of the n — 1 solution and the mass per 
winding number of the n = 2 solution A.E = E(n = 1) — E(n = 2)/2 is shown 
as function of 7 for the first (circles) and second (squares) branch arising in the 
EYMHD system in the Volkov limit a 2 = Also shown is AE for the YMHD 
solutions (triangles). Note that for the EYMHD solutions, the 7- axis represents a 
rescaled a-axis. 
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